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We have conducted a search for extended energy deposition trails left by ultra-relativistic magnetic 
monopoles interacting in Antarctic ice. The non- observation of any satisfactory candidates in the 
31 days of accumulated ANITA-II (Antarctic Impulsive Transient Apparatus) flight data results in 
an upper limit on the diffuse flux of relativistic monopoles. We obtain a 90% C.L. limit of order 
10~ 19 (cm 2 s sr) _1 for values of Lorentz factor, 7, 10 10 < 7 at the anticipated energy E to t = 10 16 
GeV. This bound is stronger than all previously published experimental limits for this kinematic 
range. 



I. MAGNETIC MONOPOLES 



The search for magnetic monopoles has been the fo- 
cus of concerted experimental effort since Maxwell's orig- 
inal formulation of electromagnetism 150 years ago[lj. 
The non-observation of an obvious partner to electric 
charge is somewhat problematic, although inflation pro- 
vides a mechanism for diluting the primordial monopole 
abundance to miniscule densities in the current epochQ. 
Searches have spanned a wide range of monopole masses 
and Lorentz gamma; all searches assume the Dirac angu- 
lar momentum quantization condition relating the unit 
of electric charge e to the unit of magnetic charge <?Mj 
via legujc = h\^. Thus far, no report of magnetic 
monopole detections $® have been substantiatedjfl-Q . 
In 1970, Parker pointed out that the abundance of mag- 
netic monopoles is constrained by the requirement that 
magnetic monopole currents be insufficient to deplete 
galactic magnetic fields @. Only in the past decade have 
astrophysical experiments been able to improve upon the 
original Parker flux bound of ~ 10~ 15 (cm 2 s sr) _1 [l(J. 
The first observational astrophysics experiment to obtain 
limits stronger than the Parker bound was MACRO. For 
monopole velocities v = (3c > 0.99c, MACRO obtained 
a flux upper limit of 1.5 x 10 _16 (cm 2 s sr) _1 . Upper 
bounds of this order of magnitude were also reported 



for 4 x 10~ 5 < [3 < 0.99[l0j]. Subsequently, stronger 
bounds were reported by the PMT-based neutrino tele- 
scopes AMANDA [ill and Baikal [H [HI, although these 
latter searches target values of 7 somewhat lower than 
for the search described herein. 

The SLIM experiment at the Chacaltaya High Alti- 
tude Laborary in the mountains of Boliva offers sensitiv- 
ity to so-called Intermediate Mass Monopoles (IMMs). 
This nuclear track detector experiment is designed es- 
pecially to search for monopoles (mass 10 5 — 10 12 GeV) 
over a wide range of velocities (including /3>4xl0~ 5 
for 1-Dirac-charge monopoles). SLIM's latest monopole 
flux limit at j8 « 1 is 6.5 x 10~ 16 (cm 2 s sr) _1 for 
Earth-crossing monopoles, or 1.3 x 10~ 15 (cm 2 s sr)" 1 
if upgoing monopoles are absorbed by the Earthfl4j. 
Stronger flux limits based on astrophysical considera- 
tions (such as an "extended Parker bound" ) of less than 
3 x 10 _22 (cm 2 s sr) _1 for IMMs, using an updated model 
of galactic magnetic fields[lj| have also been proposed. 

Relativistic IMMs 

Because of their moderate mass, IMMs may acquire 
highly relativistic velocities. Wick et al.\u^ use a model 
of monopole traversal of intergalactic magnetic fields 
(similar to the model underlying the Parker bound) to 
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estimate that IMMs created in the early universe would 
now have typical kinetic energies on the order of 10 16 
GeV. PeV-mass monopoles would therefore reach Lorentz 
factors 7 10 10 . The fact that IMMs acquire such 
"ultra-relativistic" 7 values provides a mechanism for 
their detection. Any particle traveling through a medium 
loses energy, but ultra-relativistic char ged particles do so 
dramatically, initiating bright showers [17] . As noted by 
Wick et aLpjl, experiments based on in-ice radiowave 
detection of showers are particularly well-suited to ultra- 
relativistic monopole detection because of a combina- 
tion of large effective volume and favorable scaling with 
energy. Most recently, the RICE experiment trans- 
lated their non-detection of radio emissions from com- 
pact neutrino-induced in-ice showers into a limit on the 
relativistic monopole flux using five years of accumulated 
datafli]. This was not done by a dedicated search, per 
se, but rather by calculating the efficiency, from Monte 
Carlo simulations, for monopole-induced showers to pass 
their neutrino search requirement criteria. It is through 
detection of such bright showers that ANITA is sensitive 
to magnetic monopoles. 



are then split into a 'trigger' path consisting of a heirar- 
chy of trigger conditions ("levels") and a 'digitization' 
path; if all levels of the trigger are satisfied, the digitized 
signals are then stored to disk. 

B. Summary of ANITA missions to date 

Hadronic and electromagnetic showers resulting from 
in-ice neutrino interactions produce a coherent, radio fre- 
quency Cherenkov radiation signal. Two one-month long 
missions (ANITA- 1; Dec. 2006- Jan. 2007 and ANITA- 
II; Dec. 2008- Jan. 2009) have yielded the world's-best 
limits to the ultra-high energy (UHE) neutrino flux in 
the energy range to which ANITA is sensitive [lH, [23| . 
A recent analysis of the ANITA- 1 data sample has also 
provided a statistically large (16 events) sample of self- 
triggered radio frequency signals attributed to geosyn- 
chrotron radiation associated with cosmic-ray induced 
extensive air showers (EAS) 24]. The analysis described 
herein is based on the ANITA-II data sample. 



II. ANITA 

The Antarctic Impulsive Transient Antenna (ANITA) 
is a balloon-borne antenna array primarily designed to 
detect radio wave pulses caused by neutrino collisions 
with matter, specifically ice. The basic instrument con- 
sists of a suite of 40 Seavey Corp. quad-ridged horn 
antennas, optimized over the frequency range of 200- 
1200 MHz, with separate outputs for vertically vs. hor- 
izontally incident radio frequency signals, mounted on 
a high-altitude balloon. From an elevation of ^38 km, 
the balloon scans the Antarctic continent in a circumpo- 
lar trajectory. After launching from McMurdo Station, 
ANITA-II was aloft for a period of 31 days with a typical 
instantaneous duty cycle exceeding 95%. 



A. Detector Description 

Details on the ANITA hardware, as well as the trigger- 
ing scheme crucial to the analysis described herein, are 
provided elsewhere 18[ . We here provide a brief summary 
of the hardware elements essential to this analysis. 

The front-end Seavey quad-ridge antennas provide the 
first element in the ANITA radio wave signal process- 
ing chain. These horn antennas have separate Vertical 
(VPol, defined as the zenith direction) vs. Horizontal 
(HPol) signal polarizations, and each have a field-of-view 
of approximately 60° full- width-half-maximum (FWHM) 
in both azimuth angle (</>) and zenith angle (0). Two con- 
tiguous antennas in azimuth define a logical 'phi sector'. 
Following the Seavey antennas, initial signal conditioning 
restricts the frequency bandpass to that desired for the 
primary neutrino search, namely, 200-1200 MHz. Signals 



III. MONOPOLE ENERGY LOSS IN MATTER 

Essential to this monopole search is a reliable simula- 
tion of in-ice showers caused by monopoles traversing the 
ice sheet. Our model of monopole energy loss is based on 
the muon/tau energy loss model of Dutta et aL[2||. In 
this model, energy loss by a muon traversing a medium 
is expressed as 



dE 
dx 



(1) 



The a term is the energy loss per distance due to ion- 
ization of the medium. The (3 term[26] is the sum of 
three terms reflecting bremsstrahlung, pair production, 
and photonuclear effect energy losses. The values of a 
and of the various contributions to /3 are only weakly 7 
dependent. 

Although energy loss due to ionization can be treated 
as smooth and continuous with little loss of accuracy, 
the stochastic fluctuations in pair production and pho- 
tonuclear energy losses must be explicitly modeled; such 
discrete processes, in fact, provide the catastrophic show- 
ers to which ANITA-II is most sensitive. The Dutta 
et al. model expresses these energy losses in terms of 
partial interaction cross sections with respect to inter- 
action energy, making it possible to isolate the expecta- 
tion number of particle/medium interactions at a given 
energy and replace it with a random number of interac- 
tions drawn from the appropriate Poisson distribution. 
To convert the stochastic model of muon energy loss to 
a model of magnetic monopole energy loss, the muon 
mass must first be replaced by the magnetic monopole 
mass. For monopoles, bremsstrahlung is negligible and 
is disregarded p^j. Next, due to Dirac's quantization con- 
dition, a magnetic monopole of 1 Dirac charge will lose 
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FIG. 1: Total energy loss versus 7 for 10 16 GeV monopoles, 
showing stochastic variation averaged over 50 m intervals. 
Lines show average contributions from different processes. 



energy equivalent to an electric charge of l/(2a) times 
the proton charge [17|. Accounting for this large effective 
charge only requires multiplying the expectation number 
of interactions by l/(2a) 2 w 4700. 

Figure [T] shows various contributions to the energy loss 
of a 10 16 GeV monopole. (The monopole rest mass is 
constrained to vary inversely with gamma such that total 
energy is fixed at the reference energy of 10 16 GeV.) The 
curves indicate average energy loss due to the three prin- 
cipal mechanisms, while the points show actual stochas- 
tic energy loss (as averaged over a 50 m interval in our 
simulation) . The photonuclear effect is the dominant en- 
ergy loss mechanism at 7 > 10 4 , while ionization en- 
ergy losses dominate below this value. Because the pho- 
tonuclear mechanism results in hadronic showers gener- 
ated by nuclear recoils, we may ignore Landau-Migdal- 
Pomeranchuk (LPM)[23| effects. 



IV. MONOPOLE SIMULATION 

The monopole energy loss model described above has 
been incorporated into a Monte Carlo simulation of a 
magnetic monopole traveling within the ice volume vis- 
ible to ANITA. The Monte Carlo simulation randomly 
generates a monopole trajectory (uniform over An) and 
energy loss according to the previously described model, 
then determines the voltage response of each of the 
ANITA antennas which participate in the event trigger, 
as summarized below and elsewhere. Earth curvature ef- 
fects, which are considerable at the largest radii (~600 
km), are incorporated into the model by a suitable ad- 
justment of the refracted ray geometry, relative to the 
balloon, for a given source location. 



A. Passage through Earth 

As can be inferred from Figure [TJ a down-coming 
monopole having 7 ~ 10 9 will easily traverse the entirety 
of the Antarctic ice sheet without substantial degrada- 
tion of energy, whereas a normally-incident up-coming 
monopole at those 7 values will typically range out in 
the Earth before reaching the South Pole. Upcoming 
monopoles with larger 7 values will generally reach the 
detector, although their energy loss in-transit through 
the Earth can be substantial. For calculating this energy 
loss, the terrestrial density integrated over the length of 
the monopole's path [g/cm 2 ], as a function of approach 
angle, is taken from the Preliminary Reference Earth 
Model M. 



B. Propagation Through Ice 

After calculating the energy lost by a monopole en 
route to the ice target (zero if the monopole arrives from 
above), the monopole's interaction with the ice itself is 
simulated, over discrete step sizes of 15 cm. This selec- 
tion of step size is motivated by several considerations: 
a) the time delay between signals received at the ends 
of the step should be smaller than the time scale of the 
trigger logic, which is determined by the 7 ns duration 
of the front-end tunnel-diode signal integration, b) the 
time delay between successive steps should not be larger 
than the ANITA digitization time of 0.4 ns, c) the steps 
should not be so fine that the simulation processing time 
(and array sizes) become prohibitive. For the most unfa- 
vorable geometry (a shower receding from ANITA along 
the line-of-sight), the time delay between successive steps 
is the sum of the time for the monopole to move 15 cm 
in ice plus the additional signal transit time in-air, or 
approximately 1 ns. However, this geometry results in 
a Cherenkov pattern undetectable at the payload - sim- 
ilarly, monopoles entering the ice at normal incidence 
from above (air — > ice) do not illuminate the balloon. 
Signals typically result when an up-coming monopole 
transits from rock to ice or a down-coming monopole ap- 
proaching the balloon is incident on the air-ice interface 
from above at a near-glancing angle, as illustrated in Fig- 
ureO At the Cherenkov angle of approximately 1 radian, 
the latter results in zero time delay between successive 
steps, growing as the azimuthal angle relative to the bal- 
loon increases from zero. The former results in a maxi- 
mum time delay of approximately one time bin between 
successive steps. The geometric aperture of ANITA to 
monopoles traversing the ice can be crudely estimated 
as follows. For a Cherenkov cone developing in the ice, 
the solid angle illuminating the balloon is approximately 
the product of the vertical width in 9 (sin 6 d6 ^0.2 rads 
at the lowest frequencies accessible to ANITA) times the 
horizontal width in <fi (d<f> ~0.5 rads), giving a maximum 
geometric acceptance of 0.1 /47r, or approximately 1%. 
This provides the greatest limitation on ANITA's sensi- 
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Simulated Monopole Signal 




FIG. 2: Path geometry of monopoles coming from 
above (air— Mce) or below (rock— Hce) illuminating ANITA. 
Red/pink arrows indicate monopole trajectories which can 
result in a significant signal measured at the balloon. For 
isotropic monopole flux, the fraction of monopoles satisfying 
this geometry is of order ~1%. (For the sake of illustration, 
ray curvature through the firn has been neglected, although 
this feature is incorporated into our simulation.) 



tivity to monopole signals. 



C. Simulated Antenna Response 

A Monte Carlo simulation of the radio frequency sig- 
nals caused by cascades (discussed in Rcf.|29j]) previously 
developed in the context of RICE's high-energy neutrino 
flux studies has been adapted for application to ANITA. 
In addition to the stochastic energy loss processes enu- 
merated previously, we also simulate ray-tracing through 
the firn, as well as signal loss in transmission across the 
snow surface into the air, as a function of incident angle 
on the ice-air interface. 

A typical example of a simulated monopole time- 
domain voltage profile, for a monopole simulated at 
7 = 10 9 , is shown in Figure [3J In contrast to temporally 
compact showers from neutrinos, monopoles deposit en- 
ergy over a timescale considerably longer than the min- 
imum ~400 ns time required to fill four available data 
acquisition event buffers. Schematically, the monopole 
signal can be treated as a sequence of 'subshowers', each 
having a characteristic energy deposition geometry and 
time delay relative to ANITA. As demonstrated by labo- 
ratory measurements [13, HI|, the frequency-domain radio 
wave signal due to a shower in a dense medium is typi- 
cally broad and, in contrast to typical narrow-band an- 
thropogenic backgrounds, extends over hundreds of MHz. 
The complete simulated voltage profile V(t) at each an- 
tenna is determined by coherently summing the voltage 
contributions of various subshowers in each time bin. The 
signal phase at the emission point will vary slowly with 
viewing angle (32j: however, in most cases the dominant 
emission arises from a coherent region along the track 
centered around the Cherenkov point and subtending a 
small viewing angle. Accordingly, we ignore all signal 
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FIG. 3: Simulated voltage vs time (as measured at the data 
acquisition system) in a single ANITA vertically-polarized 
(VPol) antenna channel; 7 = fO 9 . 



phases other than those from travel time in our analysis. 

Similar to showers from neutrinos, the signal strength 
is typically at least an order of magnitude stronger along 
the vertical polarization (VPol) axis compared to the 
horizontal polarization axis, given the excellent cross- 
polarization isolation (~20 dB) of the Seavey antennas. 



D. Trigger Simulation 

The heirarchical ANITA-II trigger allows a relatively 
low neutrino detection trigger threshold, while maintain- 
ing a tolerable (~Hz) thermal noise data rate written 
to disk. In contrast to ANITA-I, only signals from the 
VPol channel of the dual-polarization horn antennas con- 
tribute to the ANITA-II trigger. Following the antenna, 
signals routed through the 'trigger' (vs. 'digitization') 
path are tested for their spectral power in four frequency 
bands, approximately covering the intervals (in MHz) 
200-^350, 330->-600, 630-^1100 and 150-^1240, respec- 
tively. This partitioning is performed in order to provide 
rejection power against narrow-bandwidth anthropogenic 
backgrounds, while retaining high efficiency for sharp 
duration (temporally), large-bandwidth neutrino signals. 
The frequency-banded signals are then passed through 
a tunnel-diode, which integrates roughly 7-ns units of 
data and provides a unipolar (negative) output pulse. 
The lowest- level trigger (L0) requires signal in one of the 
four frequency bands at a level exceeding approximately 
2.3ov, with cry the rms of the typical tunnel-diode out- 
put voltage at this point. The LI trigger required two of 
the three frequency bands plus a full-band trigger within 
a 10 ns window. (In an improvement over ANITA- 1, the 
different group delays of the bandpass filters from the 
single-band stage were compensated digitally to ensure 
synchronization.) The LI triggers were then combined 
for the L2 trigger stage. For the upper and lower anten- 
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nas, an L2 was issued in one of the 16 phi sectors when 
two of the three antennas either on or adjacent to that 
phi sector triggered within approximately 10 ns of each 
other. For phi sectors containing a nadir antenna, an 
L2 was issued for each LI, and for those without, an L2 
was issued when either adjacent nadir antenna received 
an LI. Finally, the fourth (L3) stage required two of the 
three antenna rings (upper, lower, and nadir) to have an 
L2 trigger within approximately 10 ns. Satisfaction of all 
trigger tiers initiates digitization of all antenna channels 
and occurs at a rate of approximately 10 Hz. By con- 
trast, the LI and L2 rates are approximately 2-3 MHz, 
and 2.5 kHz per antenna pair, respectively. Laboratory 
studies have demonstrated that trigger efficiency differ- 
ences between ANITA- 1 and ANITA-II result primarily 
from the removal of the HPol antenna channels from the 
trigger logic, and the modified frequency banding. 
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FIG. 4: Minimum time between successive triggers (in-flight 
data). 



V. ANITA EXPERIMENTAL MONOPOLE 
DETECTION STRATEGY 

As a monopole interacts with the ice, it sheds energy 
stochastically, so the monopole signal will be comprised 
of a string of impulsive subshower events. When ANITA 
triggers, the signal is temporarily stored in one of four 
buffers and is then read to hard disk. As soon as one 
buffer is filled, it is read out and then reset for the next 
event. Since this process requires some minimum time, if 
a series of triggers occur in rapid succession (e.g., within 
1 /is), there is insufficient time to maintain at least one 
free buffer. Once all four buffers fill, read, clear, and 
reset commands are sequentially issued, requiring a min- 
imum time of tens of milliseconds, much longer than the 
monopole signal would be visible to the balloon. The 
experimental monopole signature registered in ANITA is 
therefore expected to consist of the first four threshold- 
crossings (~500 nanoseconds total) once the monopole 
comes into view; the remaining signal produced by the 
monopole ionization trail is not registered during the 
dead-time required after the four buffers initially fill. The 
500 nanosecond estimate also includes a small delay be- 
tween successive triggers, as they fill the available buffers. 
The experimental in-flight data show a minimum time 
between triggers T m j n =112 ns (Fig. [4j, to be compared 
with the 100 ns of data recorded in a standard waveform 
capture. This value is commensurate with laboratory 
measurements prior to the flight. 



VI. MONOPOLE SELECTION AND EVENT 
YIELD 

Our monopole event selection criterion therefore con- 
sists of four 'fast' triggers over a time interval T 4 ; it re- 
mains only to define how quickly the triggers must be 
registered (T max ) in order to be considered a monopole 
candidate. To minimize bias, T max was determined using 
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FIG. 5: T4 distribution for Monte Carlo monopole signal sim- 
ulations compared with DB data (points). Note that, for our 
simulations, we conservatively use T m i„ = 120 ns. Also over- 
laid is a fit to the DB data showing the extrapolation into 
the region T4 <500 ns, used to estimate the expected anthro- 
pogenic background for the signal sample. 



a 'dedicated background' (DB) subset of the total ANITA 
dataset, defined as all data taken when the balloon was 
within 300km of McMurdo base. This subset comprised 
approximately 15% of the total data. For this dataset, 
Figure [5] shows the T 4 distribution, with the monopole 
MC simulated distribution overlaid, as well as a polyno- 
mial fit to the DB distribution. For large values of 7, 
the stochastic nature of the monopole energy loss tends 
to induce threshold crossings at a rate much faster than 
that at which the instrument can trigger. Since the filling 
of all buffers will immediately incur tens of milliseconds 
of subsequent dead time, the simulated T 4 distributions 
therefore exhibit pile-up at the minimum possible value 
of 3T m i n . 

From the DB sample, we set the requirement that 
T4 < T max (= 500 ns), which eliminates all the back- 
ground event triggers. This is the main monopole selec- 
tion criteria then applied to our 'search' sample. 
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A. Trigger Efficiency 

Monte Carlo simulations of monopoles isotropically il- 
luminating the Antarctic ice sheet are used to determine 
the T4 <500 ns efficiency. Although the bulk of high- 
energy monopoles which produce four triggers also sat- 
isfy the T4 <500 ns criterion, the limited geometric ac- 
ceptance for producing four triggers (<1%, designated 
as "e a ") provides the main limitation on the total effi- 
ciency (designated as "e*"). The ratio of events satisfying 
both the simulated trigger (which includes aperture), as 
well as the T4 <500 ns requirement relative to the total 
events generated in simulation for each 7 value (typically 
20,000), are as follows: 



^0510(7) 


8 9 10 11 12 13 


e a (xKT 3 ) 
e t (xHT 3 ) 
a(e t ) (xlO" 3 ) 


0.43 3.85 4.35 7.19 7.69 
.06 1.65 3.17 4.66 6.12 
.04 .34 .48 .99 1.20 



Also shown in the last row are the errors on the total 
efficiency; our eventual limits will later be degraded by 
one statistical error bar in each energy bin. 

B. Background estimate in search sample 

After establishing the requirement that T 4 <500 ns 
based on the DB sample, the trigger time distribution of 
the remainder of the data set (the monopole 'search' sam- 
ple) was considered. It is desirable that the number of 
events passing this requirement in the search sample be 
small enough that they can be individually hand-scanned 
and their candidacy assessed against Monte Carlo simu- 
lated monopole templates. Extrapolating the DB T4 dis- 
tribution into the 'signal' region below 500 ns (Fig. [S]) 
indicates a projected yield of 0.46 events; correcting this 
by the ratio of livetime in the search sample relative to 
the DB samples (0.85/. 15) results in an a priori expec- 
tation of 2.6 background events in the search sample. 

Thermal triggers, overall, dominate the accumulated 
ANITA-II data sample and are responsible for the 'hump' 
to the right in Figure [5] Given the typical L3 trigger rate 
of 10 Hz, the probability of registering 4 uncorrelated 
thermal triggers within 500 ns can be estimated to be 
~ (5x 10 -6 ) 3 per second, implying a negligible expected 
thermal event fake yield in our search sample. 

C. T4 distributions in search sample 

The T 4 distribution for the entire ANITA-II flight is 
shown in Figure [6j We do note dissimilarities, at small 
values of T 4 , in the shapes of the DB vs. search sample T4 
distributions. This is, however, not necessarily surpris- 
ing, since our DB sample is restricted to data taken in 
the vicinity of McMurdo Station, while the search sample 
includes data taken over the entire continent. 




Time to Register Four Triggers (ns) 

FIG. 6: Time between first and fourth event trigger (T4) 
distribution for ANITA-II monopole search sample; the 
monopole event selection requirement (T4 <500 ns) is indi- 
cated by vertical line. 



In the entire data sample, only four events contain 
four rapid triggers which satisfy the 500 ns maximum 
total trigger time criterion. All four of these events 
were already classified as background, using the rejec- 
tion criteria developed for the primary ANITA-II neu- 
trino search[23j]. Specifically, the first two events were 
rejected on the basis of their temporal near-coincidence 
(occurring within 15 seconds of each other) as well as 
having source locations (after filtering for strong contin- 
uous wave [CW] components) consistent, to within ~1 
degree, with the South Pole. The third event was recog- 
nized as payload noise, on the basis of an anomalous DC 
offset, and saturated waveforms in all phi sectors. The fi- 
nal event was rejected on the basis of its power spectrum, 
which is observed to be dominated by one narrow CW fre- 
quency without which the event would be sub-threshold 
for analysis. Removing that one line from the frequency 
spectrum results in a time-domain waveform which also 
reconstructs to a known background source location in 
West Antarctica and further affirms this as a background 
event. The time-domain, as well as frequency-domain 
distributions for the four highest-amplitude channels in 
that last event are displayed in Figure [7] We note sim- 
ilarity in character between that event from the search 
sample and a typical event with T4 «500 ns from the 
DB sample (Figure [8|). In fact, the power spectra of all 
the events with T4 <500 ns, in both the DB as well as 
search samples, are found (o posteriori) to be dominated 
by CW lines. This is in marked contrast to monopole 
signals, which exhibit the broadband characteristics of 
temporally-sharp, coherent radio wavelength signals. 

Elimination of these four events in the search sam- 
ple (roughly consistent with our a priori expectation of 
2.6 background events from the DB sample) results in 
zero candidate monopole triggers in the ANITA-II data 
set. Since no satisfactory event candidates are found, 
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FIG. 7: Time domain (top four panels), as well as frequency 
domain (bottom four panels) voltage distributions for the four 
channels in one of the four search sample events passing the 
T4 <500 ns requirement. 




FIG. 8: Time domain (top four panels), as well as frequency 
domain (bottom four panels) voltage distributions for the four 
channels in one of the DB events with T4 ~500 ns. 



we therefore choose to set a limit on the monopole flux. 
However, having invoked criteria (base rejection) used 
for the primary neutrino search, we also must later de- 
grade our monopole detection efficiency by the geometric 
loss due to that anthropogenic background requirement 
(e 6 = 0.63). 



VII. CROSS-CHECK ANALYSIS 

A complete, and parallel analysis, which was intended 
to avoid the a posteriori hand scan, was performed as fol- 
lows. Given the fact that anthropogenic noise tends to be 

a) CW (continuous wave) in the frequency domain, and 

b) persisting for durations of order seconds, we initially 
restrict our candidate search sample to those events sat- 
isfying an extended CW-rejection requirement. Starting 
with the first event of a set of four filled in rapid succes- 
sion, we transform to the frequency domain and deter- 
mine the frequency bin with the highest fractional power. 
For a uniform population, and based on 256 time-domain 
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FIG. 9: T4 distribution for Monte Carlo simulations and DB 
sample, for events passing PFP < 10~ 7 requirement. Note 
the diminution of dedicated background events populating the 
T4 ~500 ns region, in comparison to Figure 5. 



samples, this quantity would be 1/(256/2), or 1/128. As- 
suming that, for man-made backgrounds, that bin will 
also have the highest spectral power for the next 3 events 
in a candidate quartet, we now calculate the product of 
the fractional power ("PFP"), for that same bin deter- 
mined from the first captured event, for all four events 
in a quartet. For thermal noise, this quantity peaks at 
around 10 -9 , for Monte Carlo simulations, this distribu- 
tion peaks at around 10 -8 , primarily because of the fact 
that: i) monopole candidates are generally close to the 
thermal floor, and ii) the ANITA trigger requires power 
in multiple trigger bands. Based on the Monte Carlo dis- 
tribution, we set a cut on this statistic at PFP < 10~ 7 , 
which corresponds to approximately 95% efficiency and 
high background rejection of the DB sample (Figure O ■ 
Having set the cut based on Monte Carlo simulations 
only, we now apply this requirement to the signal search 
sample; none of the T4 <500 ns data events survives this 
requirement, resulting (again) in zero monopole candi- 
dates. Note that application of this CW-rejection would 
actually result in a stronger limit than what we finally 
quote below as we do not incur a penalty of ef, in our 
efficiency. 



VIII. LIVE TIME 

In order to calculate a flux limit, the time ANITA 
is sensitive to monopoles must be known. Recall that 
ANITA must be both active and have four free buffers 
available to record our "cluster" of events. The buffer 
depth ("BD", or the number of available buffers to fill 
with data triggers) was not always known during flight, 
so a simulation was developed to estimate the fractional 
live time contributed by each buffer state. We know 
that when an event is triggered, a buffer becomes oc- 
cupied, and the buffer depth decrements (we use the no- 
tation "BD=4" to indicate that all buffers are empty; 
"BD=0" implies all buffers are occupied and dead time 



8 



HOLDA 
HOLDB 
HOLDC 
HOLDD 
DEAD 



J 



bd y^T^r^r^f^^ 

CLR 
TRIG 



FIG. 10: Timing diagram, illustrating the sequence of fill- 
ing and emptying event data buffers. Four rapid triggers 
("TRIG") initiate data HOLD in buffers A, B, C and D, and 
decrements the Buffer Depth (BD) variable. Once all buffers 
are filled, DEAD time is incurred during the clear (CLR) and 
reset cycles. 
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FIG. 11: Comparison of the simulated readout time model 
with data. "Simulated" quantities refer to outputs of the live 
time model, as the effective readout time is varied. At an 
effective readout time of 73 ms, the simulated livetime (blue, 
short dashed) matches the known livetime (cyan, dotted), im- 
plying a BD=4 livetime slightly less than 10 6 s. 



IX. SYSTEMATIC ERRORS 



is incurred). When an event is read to disk, a buffer is 
freed and BD increments. Figure [10] is a graphical repre- 
sentation of event acquisition and readout. 

Our simulation provides an estimate of the amount of 
time spent in each buffer state Tbd=u where i indicates 
the buffer depth; the monopole livetime therefore corre- 
sponds to Tbd=a- The only variable in the model is the 
unknown effective readout time T reac iout, which is deter- 
mined as follows. The simulation reads through the en- 
tire list of trigger times recorded in data. When an event 
is triggered, our model increments BD and, after T rea dout 
has elapsed, BD decrements. If, in the simulation, a trig- 
ger occurs when BD=0, the event is simply skipped and 
BD is left unchanged. Given the typical ANITA-II data 
trigger rate of 10 Hz, any gaps in triggers greater than 
1 second are considered dead time, and following one of 
these gaps, BD is reset to 0. Despite the individual buffer 
depth states not being known, the total time ANITA was 
live (Ttotai = T BD=1 + T B d=2 + T BD=3 + T BD=4 ) is re- 
liably tracked during flight and provides a constraint on 
our model. The 'correct' T rea dout is defined as that value 
for which the simulated T to tai converged on the actual 

Figure[TT]shows the simulation results as the trial value 
of T rea dout is varied from 0.02s to 0.14s. The horizontal 
magenta line represents the known T to taV, the intersec- 
tion of this line with the cyan arrow (simulated Ttotai) 
corresponds to the case where the live time constraint 
is satisfied. Simulation matches data when T rea dout ~73 
ms. From this value, our aggregate livetime for BD=4 is 
estimated to be 9.75 x 10 5 s. 



The primary source of uncertainty is likely due to the 
systematic error inherent in our efficiency estimate, pri- 
marily due to ice properties uncertainties and our finite 
time binning procedure (which determines the coherence 
of the time domain waveform) in simulations. This is 
assessed by running our simulation with two consider- 
ably different sets of parameters. In the first set of trials, 
the monopole is tracked over a step size of 80 cm, be- 
ginning at the first point within the ice sheet when the 
monopole registers a trigger, and extending for 8192 sim- 
ulated samples. For high-energy, upcoming monopoles 
which are capable of arriving from below after travers- 
ing the Earth, the monopole is therefore tracked over 
the warmest, and most birefringent ~1.5 km of ice start- 
ing at the bottom of the ice sheet. In the second set of 
trials (our default parameters), the monopole is tracked 
over a 15 cm step size only over the top (colder ice, with 
no birefringence) 1.5 km of ice thickness, extending over 
16384 simulated samples. The 'test' parameters result in 
approximately *~23% lower efficiency than our default pa- 
rameters, largely due to the generally-smaller measured 
signal strengths emanating from the deeper, warmer (and 
significantly more radio-absorbing) ice. This extreme 
variation should bracket the expected uncertainties due 
to our radio frequencies ice properties' parameterization, 
as well as our finite time-binning procedure for signal, in 
the limit where the step size is correctly taken to zero. 

There is some uncertainty associated with the calcula- 
tion of the monopole energy loss. To assess this, we have 
compared our simulation with the monopole energy loss 
calculated independently by the mmc code (3^. which was 
originally designed for muons and subsequently adapted 
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to model monopoles for the IceCube collaboration. Un- 
like our simulation, the mine energy loss paramctrization 
is based on several different theoretical calculations, and 
also calculates energy losses independently of our simu- 
lations. Since the kinematic regime targeted by the Ice- 
Cube analysis is approximately 5-10 orders of magnitude 
lower than for this analysis, the overlap between the sim- 
ulations is limited to gamma values of order 10 5 . For 
7 = 10 5 , we find agreement between the results of the 
two simulations to within 13.5%. 

Our overall result is relatively insensitive to the 
parametrization of thermal noise in our simulation. The 
addition of such noise is included in our simulation and 
will, on average, lead to a noticeable efficiency improve- 
ment for very large simulation samples. For our lim- 
ited simulation samples, however, this effect is mitigated 
by the stochastic nature of the monopole energy depo- 
sition, which leads to a large number of in-ice showers 
with widely varying energies, and distributed over view- 
ing angles varying by up to several degrees relative to the 
balloon. 

To account for our limited Monte Carlo statistics, we 
have included an error equal to the fractional statisti- 
cal precision on our Monte Carlo-derived detection ef- 
ficiency. To account for the uncertainty in our choice 
of T4 <500 ns as the optimal monopole event selection 
requirement, we have also folded in a systematic error 
equal to the inefficiency of our T4 <500 ns event selec- 
tion requirement (10%). Added in quadrature with both 
the energy loss uncertainty (13.5%), as well as the 23% 
systematic error determined above yields a net system- 
atic error of 28%. Our final upper limit is accordingly 
degraded by this fraction. Uncertainties in our livetime 
calculation are believed to be smaller than 5% and can 
therefore be neglected in our total systematic error. We 
note that the two dominant systematic errors above are 
independent of each other - the energy loss uncertainty is 
assessed without fully propagating the signal to the bal- 
loon, while the latter uncertainty measures the error in 
our modeling subsequent to energy deposition. Note that 
the final result is somewhat robust to changes in the over- 
all scale of dE/dx - increasing dE/dx leads to brighter 
in-ice showers, but reduces the flux of monopoles that can 
penetrate to the ice sheet from below the horizon. Con- 
versely, reducing dE / dx increases the flux arriving at the 
ice sheet, but decreases the magnitude of signal strength 
on a monopole-by-monopole basis. In general, ANITA's 
sensitivity to relativistic monopoles is primarily deter- 
mined by the detector/ice geometry. For monopole tra- 
jectories that do illuminate the balloon, we expect there 
to be tens, if not hundreds of potential showers which 
can cause event triggers. In practice, only the first four 
of these are recorded by the data acquisition system. 

Overall, we have tried to take a conservative approach 
in calculating our sensitivity, including: a) tracking the 
monopole over the only the upper half of the ice sheet 
and neglecting any signal from the lower half which re- 
sults in a conservative estimate of effective area, b) calcu- 



lating sensitivity under the assumption that systematic 
errors are uncorrelated, and therefore forfeiting any nu- 
merical advantage that would be gained by taking such 
correlations into account, c) using an estimated trigger 
efficiency which is approximately 25% less efficient than 
that employed for the primary neutrino search. 

X. PARTICLE FLUX LIMIT CALCULATION 

Particle flux is reported in units of cm s sr . To 
calculate an upper bound on the diffuse monopole flux, 
the upper limit on the number of observed particles (N) , 
live time (Tu ve ), detection area (A), and the solid angle 
from which incident particles are approaching (f2) must 
be known. In this experiment, none of the candidate 
events met the acceptance criteria, so the number of ob- 
served particles is zero, corresponding to a Poisson 90% 
confidence level upper limit on the number of observed 
particles of 2.3 (we neglect possible background contri- 
butions which, if included, would tend to strengthen our 
derived upper limit). The BD=4 live time given by the 
simulation is THt, e =975,000 s. The detection area is the 
surface area of a spherical cap of radius R=680 km and 
is therefore slightly greater than the planar projection 
area of it x R 2 (by approximately 5%). Finally, we must 
include the efficiency for events to both trigger, as well 
as to pass our T4 <500 ns (e t ) requirement, and also the 
efficiency for events to pass our base rejection require- 
ments (e&) in the expression for particle flux F (Eqn. 
we take e& to be the value used for the neutrino-search 
analysis (63%), as those search criteria are invoked dur- 
ing the final phase of the analysis in our elimination of 
the four events passing the X4 <500 ns requirement from 
monopole candidacy. 



AQTi ive (e t e b ) 

Our final results are shown in Figure Q21 We obtain 
a flux limit of order 10 _19 /(cm 2 — s — sr) at 7 = 10 10 , 
improving slowly over the next few decades in monopole 
7. For large gamma values, the results for ANITA are 
considerably lower than that of any experiment to date. 
Although the bandpasses of the RICE and ANITA ex- 
periments are comparable, the ANITA analysis strategy 
is directed at monopole detection ab initio, affording a 
factor of three improvement in efficiency relative to the 
RICE analysis, which was driven by a neutrino search 
strategy. The additional enhancement in effective collec- 
tion area (~500, including losses in collection area due to 
rejecting bases) more than compensates ANITA's factor 
of ~50 smaller livetime relative to RICE. 



XL SUMMARY 

From the non-observation of highly ionizing showers 
we have derived the monopole flux upper limits shown 
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FIG. 12: Comparison of ANITA upper limit on diffuse 
monopole flux with other results. Save for RICE, other ex- 
perimental results have been extrapolated up to our sensitive 
kinematic interval. In performing this extrapolation, the lim- 
its for 7 > 10 9 have been weakened by a factor of two, to 
account for increasing Earth opacity. 



in Fig. [T2j which are on the order of 10 19 (cm 2 s sr) 1 . 
Previously, AMANDA[ll|, Baikalpjj], and MACRO Q3 
determined monopole flux limits on the order of 
10~ 16 (cm 2 s sr) -1 for /3 greater than 0.8, 0.8, and 4 x 
10~ 5 , respectively. Although the results of this study 
cover a much narrower range of /3 values than previous 
works, it is the range that is of the greatest interest for 
IMM searches. Within much of this kinematic range 
(E = 10 16 GeV; 7 > 10 9 ), monopole flux limits from 
ANITA are stronger than the limits from any previous 
astrophysical monopole search. 
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